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field 
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The structure of the superconducting order parameter in the iron-pnictide superconductor 
BaFe2(Aso.67Po.33)2 (Tc = 31 K) with line nodes is studied by the angle-resolved thermal conductivity 
measurements in a magnetic field rotated within the basal plane. We find that the thermal conductivity displays 
distinct fourfold oscillations with minima when the field is directed at ±45° with respect to the tetragonal a 
axis. We discuss possible gap structures that can account for the data, and conclude that the observed results 
are most consistent with the closed nodal loops located at the flat parts of the electron Fermi surface with high 
Fermi velocity. 

PACS numbers: 74.20.Rp, 74.25.fc, 74.70.Xa 



The fundamental mechanism that causes high-temperature 
superconductivity is a central issue in the physics of Fe- 
pnictides. Knowledge of the superconducting (SC) gap struc- 
ture is a major step toward identifying the interactions that 
produce pamng. However, despite tremendous experimental 
and theoretical efforts, the SC gap structure remains one of the 
most controversial topics in these materials. The Fe-pnictides 
have several electron and hole bands and the interband inter- 
action is very important for superconductivity. ' Nesting be- 
tween electron and hole bands promotes antiferromagnetism 
and the s±-wave gap symmetry with a sign change between 
nested bands. ^'^ In contrast, orbital fluctuations promote an 
S++ gap without a sign change. 

Recent experiments revealed that several Fe-pnictide su- 
perconductors have line nodes in the SC gap. They in- 
clude LaFePO (T^ « 6K),5 '' heavily hole-doped KFezAsa 
(Tc « 3K)i and BaFe2(Asi_2:Pj;)2.^'^ The presence of line 
nodes in the latter material is especially intriguing because the 
transition temperature of this system in the optimally doped 
region, Tc = 31 K at a: = 0.33,'" is comparable to the 
highest Tc of hole-doped (Bai_^K^)Fe2As2 (T^ = 37 K) 
whose Fermi surface (FS) is fully gapped in the SC state i'*''^ 
In BaFe2(Asi_a;Pa;)2, substitution of isoelectronic P for As 
suppresses the magnetic order and induces superconductivity 
without charge doping. This system is very clean as demon- 
strated by the observed quantum oscillations'-' and by the 
low critical current density in the vortex state." In addition, 
the x-dependence of several physical quantities is consistent 
with a quantum critical point near the boundary of magnetism 
and superconductivity.'^ Therefore this system is suitable for 
studying the intrinsic nature of the high-Tc superconductivity 
and for determining the detailed SC gap stiTicture. 

Here we report angle-resolved measurements of the ther- 
mal conductivity k in a rotated magnetic field H, which is a 
bulk probe of the position of the gap nodes in unconventional 
superconductors. '^^^ We observed distinct fourfold modula- 
tions under field rotation within the ab plane. Combined with 
the results of the superfluid density ps,'^ specific heat C,^^ the 



band structure calculation and guided by comparison with the 
theoretical predictions for angular variations of k for different 
FS shapes and nodal structures, we conclude that the nodal 
lines form closed loops at the flatter part of the outer electron 
pocket near the intersection of the electron pockets with the 
T-X line. 

K was measured for BaFe2(Aso.67Po.33)2 single crystals — 
by the steady-state method with the thermal cuiTent q along 
the tetragonal a axis [Fig. [Ha)]. We rotate H within the ab 
plane with a misalignment of less than 0.02° from the plane. 
The temperature and field dependence of k/T of the present 
crystal agree with the data in Ref. 8, including finite residual 
K/TatT OK in zero field and ^/H dependence of k/T, 
expected of linear nodes. 

Figure[na) depicts the angular variation of where cj) 
is the angle between the a axis and H. We measured 
by changing (j> after field cooling at (p = 0°. The subse- 
quent measurement with an inverted direction of rotation did 
not produce any hysteresis in At (j> — ±45°, 

shows minima, indicating the presence of fourfold modula- 
tion. As shown by the solid lines, k{4>) at each T and H can 
be decomposed as = kq + ii2<p ± ^4^, where kq is a 

constant, and K2<j, = C24, cos 2(j) and HA(j) = C4(p cos 4(/) are 
the terms with twofold and fourfold symmetry respectively. 
The twofold contribution appears because of the difference 
in transport with the field parallel to and normal to the heat 
current, and is present even for an isotropic gap. However, in 
contrast to most previously studied nodal superconductors 
the amplitude of the twofold term is smaller than that of the 
fourfold term [Fig.[TJb)]. As we see below, this imposes con- 
straints on the shape of the superconducting gap. 

Figures[nc)-l(h) display K40/k„ at low T and low H, 
which is directly related to the SC gap structure. The normal- 
state Kn value is obtained from the residual resistivit}^ by 
using the Wiedemann-Franz law. We observe clear fourfold 
oscillations with the amplitude |C40|/K,i^O.Ol-O.O2. Such 
oscillations are due to the creation of unpaired quasiparticles 
in the regions with small or vanishing gaps by the applied 
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FIG. 1: (Color online) (a) Angular variation of «:(<^) in 
BaFe2(Aso.67Po.33)2 with a fit to = kq + + K4</, (blue 

line). The green (yellow) line represents the twofold (fourfold) term 
obtained by the fitting. The inset illustrates the measurement config- 
uration, (b) The magnitudes of K24, and /t40 at 2 K are plotted against 
the field magnitude H. (c)-(e) The fourfold term K4,^ at fioH = 2T. 
(f)-(h) at T = 2 K. 

field. Previous k/T{T, H) and Ps{T) measurements^ already 
indicate the presence of nodes rather than the gap minima, so 
we focus solely on the location of the nodes. The strength 
of the pair breaking in a given near-nodal region vanishes if 
the Fermi velocity vp ts parallel to H, and is maximal for 
Di? _L H . At the lowest energies this can be explained by 
the Doppler shift of the quasiparticle energy by the amount 
Se{r,VF) = VF{k) ■ Ps{r), where Ps is the superflow due to 
the vortices in the plane normal to i?,"' and more microscopic 
theories retain this structured Hence rotating H provides a 
probe of the node positions. 

BaFe2(Aso.67Po.33)2 is a multiband system. We obtained 
the FS from the band structure calculations including spin- 
orbit coupling using a density functional theory implemented 
in the WIEN2k code. The FS consists of five quasi-cylindrical 
pockets, three hole pockets at the center of the Brillouin zone, 
and two electron pockets centered at its corners [Fig.|2la)]. 
Because of the body-centered tetragonal symmetry, the 
cuts of electron sheets are elliptical, with only twofold sym- 
metry. These ellipses are rotated by 90° between kz=0 and 

±— , forming a so-called "snake that swallowed a chain" 
I— i- 

shape. As shown in Figs.[3la)-3(c), flat regions appear in the 
electron pockets around point A (along T-X) and the FS has 
a higher curvature near point B. 




FIG. 2: (Color online) (a) Fermi surface of BaFe2(Aso.67Po.33)2. 
The shading represents the in-plane Fermi velocity v"}^ . Note that 
the flat parts of the outer electron sheets [near point A in Figs. 3(a)- 
3(c)] have high v^^ values, which correspond to the gray shading 
in Figs. 2(b) and 2(c). (b), (c) The DOS, N{£f) (oc C/T) and 
N{eF){vf'Y (« k/T) of each band as a function of v'^ . 



The first question is at which FS(s) are the nodes located. 
We believe that the nodes are on the electron sheets. First, re- 
cent bulk-sensitive laser angle-resolved photoemission spec- 
troscopy (ARPES) experiments report isotropic gaps in all 
three hole pockets around Z point. This by itself excludes 
the rf-wave gap symmetry and we focus on the Aig represen- 
tation. Second, there are significant changes in the shape of 
the hole FSs with x in BaFe2(Asi_a;Pj;)2, while the shape of 
the electron surfaces remains nearly the same.^" The slope of 
the T-linear penetration depth, indicative of robust line nodes, 
simply scales with Tc in a wide range of xP implying little 
change mvF at the near-nodal regions. Finally, in an applied 
field the residual k/T exhibits a clear ^/H dependence^ over a 
much wider range than the residual specific heat C/T, which 
looks nearly linear above 5 T.— This suggests that the nodes 
are located in regions with higher in-plane Fermi velocity v'^p, 
as those give a lower relative contribution to C/T [propor- 
tional to the density of states, DOS, N{eF, H)] than to trans- 
port [proportional to N{eF)iv'j:^)^]. The shading of Fig.|2a) 
represents the magnitude of v'j:^ with the red (blue) area denot- 
ing regions with high (low) Fermi velocity. It is clear that the 
hole pockets have low velocity and heavy mass quasiparticles, 
while the electron FS sheets have high velocity and light mass. 
Considering that the electron pockets have a lower scattering 
rate than the hole pockets, their contribution to the total ther- 
mal conductivity is even more pronounced. Note that these 
regions with high uf?" also have a higher Doppler shift Se for 
H _L ab. Together, these observations strongly support the 
scenario with the nodes on the electron sheet. 

The above analysis also suggests that the nodes may be lo- 
cated in the flat parts in the outer electron sheets, where the 
Fermi velocity is highest. Indeed, Figs.|2b), and 2(c) depict 
N{eF) (oc C/T) and N{eF){vf f (cx k/T) as a function 
of til?'. The gray shaded areas indicate the high velocity re- 
gions with v'j:^ > 0.3 X 10^ m/s, flat parts in the outer electron 
sheets. Although their contribution is less than 5% in the total 
DOS, it exceeds 30% in the net k/T, or more if the longer 
mean free path is accounted for. 

The profile of under a rotated field allows us to fur- 
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FIG. 3: (Color online) (a)-(c) Possible node positions on the electron 
Fermi surfaces. For fc^ = and ±2f , electron surface cross-sections 
are approximately elliptical, (d)-(g). (d), (e) When the nodes are on 
the flat part, the field-induced states near all nodes (filled red circles) 
have a comparable contribution to the heat current for iJ || q (d) and 
for iJ _L q (e), so that the twofold component of k is nearly absent, 
(f), (g) When the nodes are on the sharp ends, the field-induced DOS 
is small for nodes with vf nearly parallel to H (black crosses). The 
nodes with more quasiparticles (red circles) have a different angle be- 
tween Vf and q for the two field orientations, leading to a dominant 
two-fold anisotropy. (h) Comparison of for the two cases with 
the scattering rate 7 — T /2nTc = 0.01, phase shift of scattering 
5 = 60°, and the ratio of interband to intraband scattering 5v — 0.9. 



ther explore the nodal locations. The fourfold oscillations 
of of similar amplitude to that in materials with verti- 
cal line nodes argue against nearly horizontal line nodes sug- 
gested by several theories. ^^-^^ Depending on the temperature 
and field range either minima or maxima of can indi- 
cate nodal directions."^ The inversion of the oscillations re- 
ported in Fe(Se,Te) with a similar FS structure^' occurs at 
a much higher temperature and field than our range (0.03 < 
T/Tc <0.13 and 0.007< H/Hc2 <0.06). Therefore the min- 
ima at (p = ±45° observed in our low-T and low-H range 
indicate that the nodes are located at the position of the FS 
where vp is nearly parallel to [±1, ±1, 0] directions. 

In the Aig symmetry we consider three possible nodal 
structures shown in Figs.|3la)-3(c); (a) eight vertical line 
nodes, (b) closed loop line nodes in the flat part^^ and (c) 
nodal loops in the high curvature part. To distinguish be- 
tween those possibilities we compute the profile of the ther- 
mal conductivity as a function of the field angle using the mi- 
croscopic approach of Ref. 18 generalized to the multiband 
system in Ref. 22, which has been successful in describing 
C/T and k/T oscillations on equal footing. We take one elec- 
tron and one hole FS, assume the hole FS to be fully gapped. 



•^0.6 



0.4 - 



0.2 - 



0.8 - 
y'0.6 
0.4 



• 1 ' 


1 ' 1 ' 1 ' 1 

— 7=0.02 5v=().9 




— ld>.0Sh\=0.9- 


- 


leO.OS 5v=0.5 - 


(a) 




-(b) 


1 1 1 \ 


, 1 , 






-(c) TIT -0.02 
H/H ' -0.01 



0.2 0.4 0.6 0.8 1 
T/T „ 




45 
(|) (deg) 



FIG. 4: (Color online) Results of the calculations for the 3D FSs 
with nodal loops, (a) The superfluid density remains linear down to 
low temperatures for different values of the scattering rate 7. (b) At 
the same time a large residual term in the thermal conductivity in zero 
field is observed, (c), (d) Anisotropy of thermal conductivity at low 
T for several values of H. A large fourfold term with a minimum at 
45° is observed for cleaner material with sizable inter-band scattering 
(as well as for the dirtier case with small inter-band scattering). 



and consider a nodal order parameter on the electron sheet. 
An intuitive physical analysis connects the smallness of the 
twofold term with the nodes in the flat part of the ellipse — 
see Figs.[3jd)-3(g). The two ellipsoids represent the electron 
sheets at ^) and (f , — f ) points respectively. Recall that 
quasiparticles are predominantly generated at locations where 
Vp J- H (filled red circles), and quasiparticles with Vp \\ q 
contribute more to the heat transport. If the nodes are in the 
high curvature region [in the case of Fig.Qc)], the Fermi ve- 
locities at the nodes are not parallel. For the field parallel (nor- 
mal) to the heat current, the majority of unpaired states have 
Vp nearly normal (parallel) to the heat cuiTent [Figs.|3lf),(g)] 
yielding a dominant twofold component [see the red curve in 
Fig.[3jh)]. In contrast, the nodes on the flat part of the el- 
lipse [in the case of Fig.[3jb)] have a nearly identical direc- 
tion of the Fermi velocity [Figs.Qd), and 3(e)], and contribute 
equally to the q for both H\\q and H J- q [see the black curve 
in Fig.lJh)]. Thus our observation that a large fourfold com- 
ponent dominates over the twofold one showing clear minima 
near ±45° eliminates the possibility of the nodal loops at the 
sharp edges of the FS (Fig.[3jc)). It also argues against vertical 
line nodes [Fig.[3ja)] as this order parameter structure would 
essentially average the two curves in Fig.|3lh), leaving a large 
twofold anisotropy2^ not seen in our experiment. 

To confirm that the nodal loops as sketched in Fig.|3lb) con- 
sistently describe the bulk of the experimental measurements, 
we performed calculations for a three-dimensional (3D) FS 
given by 2mEf = {kxa)^/{l - ^cos^) + {kya)^/{l + 
/icos-^). For /i = 0.5 the FS closely mimics that ob- 
tained in the band structure calculations. We took the order 
parameter in the form A = AeY{(j), kz), where Y{(j), kz) = 
1 - \r{kz)\ + r{kz) cos 2(1), and r(fc^) = -0.6cos%£. We 
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also took the hole FS to be cylindrical with an isotropic 
gap A/j, and self-consistently determined both Ae(T, H) and 
Ah{T,H) assuming a dominant interband pairing interac- 
tion. Figure|4] shows that within this framework we can ex- 
plain all the salient experimental features: linear-in-T pene- 
tration depth, a large residual term in the thermal conductiv- 
ity^ and the dominant fourfold profile of under rotated 
field with a minimum along the 45° direction. We therefore 
believe that this is the most likely location of the gap nodes in 

BaFe2(Aso.67Po.33)2- 

Finally we discuss why such a nodal loop structure is re- 
alized. Experimentally, recent neutron scattering measure- 
ments suggest that the line nodes should create only a lim- 
ited area of sign-reversal on a single FS,^^ which is consistent 
with the present nodal loop structure. The Raman scattering 
in Ba(Fe,Co)2As2 suggested a similar loop structure of gap 
minima.— A number of model calculations predict that the su- 
perconducting gap in the electron FSs can be more anisotropic 
than that in hole bands, and that nodes can appear in some pa- 
rameter range. In that case placing the nodes in the flat 
FS regions with low DOS does not significantly reduce the SC 
condensation energy. In addition, it has been suggested that 
the nesting between the FS regions with the same orbital char- 



acter is important for the pairing interaction. We note that the 
change of the orbital character from xy to xz + yz occurs near 
the rim of the flat part of the electron pocket^ implying that 
the orbital character may favor the closed-loop nodes. Thus 
the pairing interaction and low DOS stabilize the nodal loop 
structure with no serious reduction of Tc. A more detailed 
fully microscopic calculation is required to clarify the origin 
of the closed nodal loop structure, which is unique among un- 
conventional superconductors. 

In summary, we determined the nodal gap structure of 
BaFe2(Aso.67Po.33)2 by angle-resolved thermal conductivity 
measurements. We found distinct fourfold oscillation in K{(f)) 
with minima at [±1,±1,0] directions. The observed results 
are most consistent with the closed nodal loops at the flat part 
of the electron FS. 
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